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ABSTRACT 

We consider systems or items which are placed on test at time zero, func- 
tion for a period, and die at some random time. Failure may be due to one of 
several causes or modes. A model for this situation is that at birth nature 
chooses a life time Y^ from a population of times until death from mode i. The 
time at which the item fails is min(Yj) and is the only life actually observed. The 
parameters of the life distributions may depend upon the levels of various stress 
variables the item is subject to. Maximum likelihood estimation methods are 
discussed in general. Specific methods are discussed for the smallest extreme- 
value distributions of life. Monte- Carlo results indicate the methods to be 
promising. Under appropriate conditions, the location parameters are nearly 
unbiased, the scale parameter is slightly biased, and the asymptotic covariances 
are rapidly approached. 

1.0 INTRODUCTION 

Suppose an item or system (a person, battery, computers, etc.) is brought 
to life at time zero, functions for a period of time, and then fails. We also 
suppose there are a finite number of possible causes (or modes) of failure 
labeled 1, 2,. . ., M. For example, people die of cancer, accident, heart 
attack, etc. A simplified conceptual model for actual item lifetime is that at 
birth nature chooses a failure time Y m from a population of failure times due 
to mode m. (The Y m may or may not be independent.) The observed life of 
the item is then min(Y m ) and it is known only that all the other lifetimes exceeded 
this value. Such a process is a form of progressive censoring. Cox[l) appears 
to have been the fir': v .o formulate this competing risk failure time model and 
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applied it specifically in the Instance of independent exponentially distributed 
failure times for two possible modes of failure. 

Herman and Patell [4] extended the model to several independent competing 
risks, developed the maximum likelihood estimators for exponential and Weibull 
distributions of failure, and derived expressions for the asymptotic covariance 
matrices of the estimators. McCool |9j considers certain confidence interval 
techniques in this situation when life has a Weibull distribution. Moeschberger flQ] 
extended the model to the case of dependent or correlated failure times. Other 
closely related developments are by Sampford [13] , Moeschberger and David [11J 
and H oel [5j. 

For any of the failure modes it may very well be that the population of failure 
times from which the item life is drawn may depend upon the level of various 
stress variables the item is subjected to during its life. For example, if we were 
to assume human life until death by cancer to follow an extreme value distribution 
with location parameter ft and scale parameter o, then n~ n( x y • • » * n ) ma y 
be a function of exposure to pollution (Xj), amount smoked per day (x^), etc. 

McCool [8] considers estimation, confidence interval, and multiple comparison 
methods when the life distribution is Weibull, failure modes act independently, 
and n observations are obtained from each population. (He essentially provides 
an ANOVA generalization.) Nelson (12] sketched a general approach to maximum 
likelihood estimation in the extension of the multiple regression situation. He 
applies the method in the case of a single stress variable (temperature), indepen- 
dent and lognormally distributed lifetimes and a number of failures by each mode 
for each level of the stress variable. 

In this report we first present the general model for the competing failure 
modes assuming the location parameters for each mode are expressible as linear 
functions of the stress variables and the failure modes act independently. We 
then present the general form of the likelihood function and the likelihood equations 
are derived for the extreme value distributions. Solving these equations using 
nonlinear least squares techniques provides an estimate of the asymptotic co- 
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variance matrix of the estimators, Jaoveral Monte- Carlo experiments were per- 
formed. The results of these indicate that, under appropriate conditions, the 
estimators approach unbiasedness and their asymptotic covariances at reasonably 
small sample sizes. 

2.0 THE GENERAL MODEL 

Assume that an experiment is performed in which I items are life tested 
at varying combinations of J stress variables and that the design is specified 
by the design matrix 




(2-1) 


The response observed for each item is a lifetime y^ (i - 1,1) and the mode by 
which the item failed, m^. The number of modes by which the item might fail 
is denoted M. 

We assume that for observation i and mode m the cumulative distribution 
function of time until failure is defined by 


F< m) |y; Ml m) , a (m ] 

and the density function by 

« l (m, [r. *! m) . ^ 

where is a location parameter and is a scale parameter. We thus 

assume that <r^ is constant for all observations. We also assume that the 
location parameter is a function of the stress variables, in particular 




ll 



( 2 - 2 ) 


In general we will not necessarily have 

(m) _ >') 


= o' 


for m 4- m' 


(2 r3) 



nor 
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S}“> - ^< m,) for m + m' 


(2-4) 


That Is, we do not a priori assume any of the failure modes to have parameter 
values in common. 

We further assume that lifetimes Y| m) are drawn at random and indepen- 
dently from Fj m ^ and that the observed lifetime is the smallest of these, min 
(y t ) . This latter assumption implies an item can fall by at most one mode. 

Nelson [12] indicates a situation in which Items may fail by one mod$, be repaired, 
and then fail again by another mode. It is also possible that testing of an item 
may be terminated before it fails. These situations may be easily included in 
the model but will not be considered in this report. 

Clearly, the profusion of arguments, subscripts and superscripts will be 
quite unwieldly. We will thus delete them whenever context clearly indicates 
(he appropriate values. 


3.0 THE LIKELIHOOD FUNCTION 

We first consider the likelihood function corresponding to the i*® 1 observa- 
tion. If that observation is a failure by mode m, then the likelihood is 



and 


1 - 





(3-3) 


That is, we multiply the density function corresponding to the mode of failure 
by the probability the item did not fail by any of the other modes. The likelihood 
is a simple product due to the independence assumption. If testing is terminated 
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without an observed failure, the likelihood is given by 

M 

*,-T 

1 



(3t4) 


The overall likelihood Is then the product of the likelihood for each observation 
and given by 

L “TTfj (3-5) 

Taking natural logs of eq. (3-5) 



(3*6) 


Is obtained where 


In L s In L ^ fif* , /3< 2: , 


a(M) 
P 1 - 


,(M) (1) (M) 

f U » * • • » V 


; y 1 yj 


(3-7) 


To obtain maximum likelihood estimators, the standard method (assuming It 
works) Is to solve the system of JM (nonlinear) equations 


9 In L _ 0 

8S{"' ) 


(m= 1, M) 


L = 0 (m= 1, M) ) 

[ 


(3-8) 


9 In L _ Q 
9o (m) 


(m = 1, M) j 


If the different failure modes have no parameters In common, as will generally 
be true by eqs. (2-3) and (2r4), this system of JM equations splits Into M separate 
systems of J equations which may be Independently solved. In particular, the 
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> (3-9) 
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Z_y a o< m > 


mode m 
failures 
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4.0 WEIBULL AND EXTREME- VALUE (SMALLEST) DISTRIBUTIONS 
Two commonly assumed life distributions are the Weibull and the smallest- 
extreme value distributions. 

The cumulative distribution function for the Weibull distribution is given by 

F(t ;Q!,/3) = 1 - exp[-(t/a)0) (4-1) 


where 

a "scale" parameter 

0 "shape" parameter 

It is easily verified that if T follows this Weibull distribution, then the 
application of the transformation Y = In T yields the smallest extreme value 
distribution with cumulative distribution function 

F[y\n,(r) * 1 - exp(-exp( (y - *t)/oj) 


with 
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a * 1/0 


M - in (Q() (4-2) 

It Is thus sufficient to consider this extreme value distribution. In particular, if 
we apply this distribution to the model previously developed we assume that for 
mode m failure and the i** 1 observation the density function for life Is 


r i v - u (m) 

fj m * y* m« * — — ex P \ — exp 

L J „(m) I „(m) 


y- 


(m) 


,(m) 


and the distribution function is 


i T 




= 1 - exp exp 


y - 


(m) 


,(m) 


(4.3) 


(4-4) 


where 


(tn) (m) +8 (m) 

Ml P* z u + • • • + Pj z ij 


(4-5) 


Hence, we obtain the partial derivatives required for the likelihood equations as 
(deleting m superscript) 

8 In f. ^ 


t , 2 V 1 


8 a a 2 


8 In f. -z 


i=^il + !iiexp y -lil 


80 


j 


8 In (1- F.) y- m ft- V 

- • ■ - * expj - 

da 2 \ a 


80i 


exp 

a V a 


(4-6) 


Substituting these into eqs. (3-9) we obtain the likelihood equations for any 


specified mode as 
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all 
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These nonlinear equations have no closed- form solution and need to be solved by 
some iterative method. 

It Is well-known that, under appropriate regularity conditions, maximum 
likelihood estimators (MLE) are consistent and asymptotically normally distributed. 
It is also known for the Weibull and extreme value distributions that the minimal 
sufficient statistics are the trivial ones consisting of the order statistics. It thus 
becomes of considerable interest to determine how rapidly the MLE's approach 
their asymptotic unbiasedness and covariance structure. The asymptotic co- 
variance structure and the Cramer- Rao lower bound are defined by the inverse 
of the Fisher information matrix (Zacks (14, p. 194] and Kendall and Stuart 
(7, p. 28]). . 

For the model considered in this report the Fisher information matrix is 
the block- diagonal matrix 


[i« 0 

0 I< 2 > 


0 


0 



(4-8) 


where 
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(4-9) 


In the following section It will be shown how an approximation to this expectation 
occurs naturally when solving the likelihood equations via nonlinear least squares. 

5. 0 SOLUTION BY LEAST SQUARES 
Using the general nonlinear regression model 

W=h(£,0)+e (5-1) 

T T 

where | * (1^, . . . , £ j) is a vector of independent variables, « = (* r • ■ , 

T 

€ n ) is a vector of random errors, 0 » (0p . . . , 0 ) is a vector of parameters 
to estimate, and the maxtrix of observed data is 


W, 


11 


; 1J 


(5-2) 


W_ 






n '’nl ' • • ‘‘nJj 

Approximating the function h as its truncated Taylor series gives (Draper 
and Smith [2, Ch. .10) ) 


h «,0) = h«,0 o ) + 



<V 


a h(g , 0 Q ) 

~ ~ 


(5-3) 
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Also, let 

0 Q = prior guess for 0 
y - correction vector 



(5-4) 


and X be the matrix of partial derivatives of h evaluated at 0 « 0 . That is, 


X u, k* 


3 h(€ u ,f>) 
a o 


0=0 


(5-5) 


With this notation eq. (5-1) can be approximated as 

W - h° * Xy ft 

which has (least squares) solution 


y - - 0 Q ) - (X T Xf 1 X T (W - h°) 


(5-6) 


(5-7) 


Thus y is the estimated correction vector. This method may be applied to the 
solution of the likelihood equations if the W's, 0's, and h of eq. (5-1) are 
defined as 


w l = 0; 0 i m Pv M^,0) 

0 /ij 


Wj.O; Uj - yfjt 

V‘ 0i ‘W" h(£ j*i, 0 »'- iL 77 


(5-8) 


With this notation we see that the partial derivatives of h (1,0) required in 
eqs. (5-5) and (5-5) are actually the second partial derivatives of the log- 
likelihood function required in the expectation of eq. (4- 9) . Thus when the X 
matrix is evaluated one has an estimate of the covariance matrix of the param- 
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eters. When used together with Monte- Carlo simulation experiments, substi- 
tution of the known parameters into the second derivatives and the observed 
sample lifetimes for each trial may be done. These results are then averaged 
to obtain an estimate of the expectation of eq. (4-9). The application of least 
squares methods to Maximum Likelihood estimation has been treated by Jennrich 
and Moore [6] . 

6.0 SIMULATION EXPERIMENTS 

Three series of sampling experiments were performed. In each of these ex- 
periments only two failure modes were considered. The first experiment in- 
vestigates the effects of censoring and sample size when the two modes have con- 
stant location parameters and identical scale parameters. The second experi- 
ment investigates the effects of various experiment designs in the case where 
mode 1 location parameter ? ■> a linear function of a single stress parameter and 
mode 2 location parameter is constant. Both modes had identical scale parameters. 
The third experiment investigates the effects of ill-conditioning on the estimators. 

For each experiment the true values of the parameters were known. Uni- 
formly distributed random variables were generated and observations for each 
mode were obtained by inverting the cumulative distribution functions. For each 
observation one lifetime for each mode was generated and the observed lifetime 
chosen as the smaller of the two. Only estimates of mode 1 parameters were ob- 
tained and recorded. The starting values for the Iterative solutions were the true 
values used to generate the observations. 

For each simulation, several pieces of information were obtained. The 
means, mean squared errors, and covariance matrix of the estimates were re- 
corded. The average amount of censoring by mode 2 was recorded along with a 
separate mean value of the estimates of the scale parameter for each amount 
of censoring. Also the Fisher information matrix of eq. (4-9) was evaluated for 
each sample and its mean value recorded. This was then inverted to obtain an 
estimate of the asymptotic covariance matrix of the parameter estimates. 

The goal of these experiments was to obtain information concerning the 
means and variances of die parameter estimates and to investigate the approach 
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to asymptotic behavior. Due to the iterative nature of the solutions required 
some sample sizes were rather small in order to avoid excessive computer times. 

The following three sections describe each of the experiments and the results 
obtained. 

6. 1 EXPERIMENT A 

This series of simulation! Investigates the performance of the MLE's for 
two populations with Identical scale parameters and location parameters of varying 
separation. In particular the parameters chosen are 

p (1) = 2.2 

H {2) = p (1) + A, A = (-0.1, 0, 0. 1, 0. 3) 

(These are highly skewed and distinctly non-normal distributions. Figure 1 shows 
one such extreme value density function and the normal density with the corres- 
ponding moments.) 

Sample sizes of n * 5, 10, and 20 were considered. Not all combinations 
of A and n were investigated. Those combinations investigated are given in 
table I. There were 1000 simulations performed for each case. Any sample 
which had less than two mode 1 failures was rejected and another sample chosen 
in its place. 

The results of these experiments are presented in table I. The first two 
columns specify the combination of 1 and n. Column three indicates the pa- 
rameter and column four the moan of the 1000 estimates for that simulation. 
Column five presents the estimated standard deviation of the mean reported in 
column four and Is obtained by dividing the observed standard deviation of that 
parameter by the square root of the number of simulations. 

The major points are that the location parameter is slightly biased toward 
lower values and that the scale parameter can be strongly biased. The bias in 
the location parameter is worst for the smallest sample sizes and the most 
censoring. The bias in the scale parameter is close to -l/n where n is the 
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sample size. The bias In <r could also depend upon A but this is not perfectly 
clear. Figure 2 plots the means of a and of p as a function of A for the three 
sample sizes. Also Included is a bar indicating plus and minus one standard 
error of the mean. The bias may decrease as A decreases. 

Harter and Moore [3J report on similar experlihents. . They were not con- 
cerned with competing failure modes but ra*her with the case where the rj 
smallest and/ or r 2 largest observations are censored. For this case they re- 
port rapid approach of the estimators to their asymptotic properties (if censoring 
is not severe) and that the scale parameter estimate is biased downward by a 
factor of l/m where m = n - r^ - r 2 is the number of uncensored observations. 
Similar results were found in the competing failure mode type of censoring. 

The last four columns of table I provide the observed covariance matrices 

a /I j /IV 

of the estimates (V ' and o' 7 and the corresponding estimate based upon asymp- 
totic theory. 

The observed covariance matrix was computed as 

v|>g 

cov|?» ^ ^ - M(^j - M<cr)j 

where 

k = 1000 trials 
and 






(6-2) 
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The estimate of the asymptotic covariance matrix was obtained as follows. Sup- 
pose the sample size Is n where y^, . . y m are mode 1 failures and y +1 , 
. . . , y n are mode 2 failures. Then 


m 




where the ju and <j are the known values of and used in generating 

the observations. It is thus possible to estimate 


2 

In L 

2 ^ 

d 6 In L 

d » 

d nd a 

9 2 In L 

9 2 In L 

^9 n 3 v 

9 <7 2 


by obtaining the average of each expression in 6-3 over all the simulated samples. 
It should be noted that for small A and n this approximation may be biased to an 
unknown degree due to rejection of samples with insufficient mode 1 failures. 

The covariances are estimated by calculating the Inverse of this matrix. 

Examination of table I shows that for A = 0. 3 the asymptotic and observed 
covariances are quite close for all n's. For A = 0. 1 the diagonal elements are 
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all quite close while the covariance is somewnat different for n = 5. Kven for 
A ■ 0 there is good agreement except for the n ~ 5 case. For A * -0. 1 it ap- 
pears the asymptotically based approximation is a considerable underestimate. 

The general conclusion is that if A a 0 the estimators approach the asymp- 
totic covariance structure very rapidly. Large A Implies little censoring by 
mode 2. These results are in qualitative agreement with Harter and Moore [3] . 

6.2 EXPERIMENT B 

This series of simulations is the simplest nontrivial situation in which the 
location parameter of one of the modes is dependent upon a stress variable. The 
models chosen for simulation are graphed in figure 3 and defined by 

p (1) = 2.2 +0.2 z 2 
p (2) = 2.5 
c^W^O.l 

The goal of the experiment was to investigate further the bias in the estimates 
and the effect of different experiment designs with respect to z 2 . In particular, 
three different sample sizes were used (n = 7, 12, and 22). For each sample 
size, three distributions of these observations with respect to z 2 were used. 

For each, z 2 was restricted to the values 0, ±1. One design concentrated all 
the replication at z g = 0, the second design spread the replication as evenly 
as possible among the three points, and the third design spread the replication 
evenly between the two extreme values. The results are based on 100 simulations 
of each case. 

Table II presents the means and the standard deviations of the means of 
each parameter. For the n = 7 sample size it appears that /3^ is biased 
slightly toward lower values, * 8 apparently unbiased, and is biased 

to lower values. The amount of bias of a does not appear to depend upon the 

A 

design. The standard deviation of 0 2 decreases as the design places more 
replication at the extremes of the range. 
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Table III compares the observed covariance matrices, the estimated asymptotic 
covariances, and the covariance matrix of the location parameters assuming naive 
normal theory. The observed and estimated asymptotic matrices were calculated 
as In experiment A. The normal theory matrix is obtained by assuming that the 
estimation process is equivalent to multiple linear regression with normally dis- 
tributed error term and ignoring the fact of competing failure modes. That is, 
the matrix is given by 

(1.645) o 2 (ZV 1 
2 

where Z is the design matrix and 1 . 645 a is the variance of the smallest ex- 
treme value distribution with scale parameter a. This describes the covariance 
structure of the location parameters only. Because of the limited number of sim- 
ulations, these results are only indicative of the true behavior and only cursory 
mention of these results is made here. 

For all cases, the observed and asymptotic covariances between and /3g 
are near zero as in the normal situation. The covariances between a and both 
/3j and j §2 are also small as in the normal case. Inspection shows the normal 
approximation to be somewhat too large in general. The asymptotic approximation 
appears to be closer for v$^) but tending to under estimate v^ 2 )' In order to 
provide more accurate comparisons, considerably more simulations should be 
performed. The results obtained here indicate that the approach to asymptotic 
covariance structure may be reasonably rapid. 

6.3 EXPERIMENT C 

This series of experiments was intended to examine what might be termed 
estlmability conditions. This is most easily described with reference to figure 4 
where p' 7 = ^ + 0 2 z + ft* are plotted. P'ive different combinations of the 
/3's were used and are presented in table IV. and illustrated in figure 4. The scale 
parameters in all cases are 

oM-oTO-0.1 

For each simulation the design consisted of n observations each at z = 

(- 2 , - 1 , 0 , 1 , 2 ). 
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With model A, there will be very few mode 2 failures at z = ±2 and it will be 
relatively clear that a second order polynominal is required to fit the data. In a 
more extreme situation such as with model E, however, there will be few if any 
mode 1 failures at z = 1, z = 2. It will then be considerably more difficult to 
determine if the model should be as in E or as in F. In essence, the censoring 
caused by mode 2 failures makes estimation of all three j3 parameters difficult 
much as if in the ordinary linear regression situation observations were not made 
at z = 1, z = 2. That is, as if we were faced with ill-conditioned normal equations. 

Tables IV and V(a) to (d) present the results of 100 simulations for various 
combinations of model and sample size. Table IV presents the means of the pa- 
rameter estimates while table V(a) to (d) presents the observed covariance 
matrices of the estimates, the estimated asymptotic covariance matrices and naive 
normal theory approximation. 

Considering first the means of the parameter estimates it is seen that cr^ 
is biased low for each sample size and the amount of bias appears to depend only 
upon the sample size, not upon the amount of censoring. 

The locus of ths location parameter estimates is more biased as there is 
more censoring by mode 2 failures as is evident in figure 5(a) to (e) . In these 
figures a solid line indicates the true value of p ^ . The graphs of the mean 
values of are indicated with dashed lines. For model A where there is 
minimal censoring, there is negligible bias. For model C where significant cen- 
soring begins to occur there is a minor but noticeable upward bias where the cen- 
soring occurs and similar bias downward where negligible censoring occurs. For 
model E where almost all of the failures for positive z are mode 2 failures, 
there is considerable upward bias at the right end of the function but minimal bias 
at the left end. These results are quite reasonable of course. Where failures 
are not observed, the maximum likelihood procedure pushes the location param- 
eter toward as large a value as consistent with the data where there are failures 
observed. 

For each simulation series, the observed covariances, asymptotic co- 
variances, and normal theory covariances of the scale parameters are given in 
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table V(a) to (d) and certain plots of the results given in figures 6 to 8. 

Table V(a) presents the results for n = 1 and cases A and B only since it 
was difficult to get batch simulation runs to complete without numerical difficulties 
in the other cases. Tables V(b) to (d) present all the results for n = 2, 4, and 8. 

A A 

The naive normal theory approximation would predict that /?g and o’ are mutually 
independent pnd independent of 0 1 and /3g. Examination of the tables indicates 
that for case A (with minimal censoring) the corresponding observed and asymp- 
totic covariances are indeed relatively small. The covariances become larger 
as the degree of censoring for positive z values increases. Normal theory 
and intuition also indicate significant covariance between and /3g. This is 
also borne out qualitatively. 

Since one picture is worth a thousand words, plots of the variances of 
are indicated in figure 6, variances of $ 2 in figure 7 and variances of $ 3 in 
figure 8. In each, the horizontal scale is n, the number of replications of the 
design, and the vertical scale is the variance of the estimator. These figures 
indicate that when censoring is infrequent (model A) both the normal theory and 
asymptotic variances behave qualitatively and quantitatively similar to the ob- 
served variances for even the smallest sample sizes. As the censoring becomes 
more severe the observed variances are much higher than either normal theory 
or asymptotic variances for small n while they become comparable for larger n. 

It is easy to understand why the normal theory variance is an underestimate. 

It is due to the (effective) inestlmability for many samples because many observa- 
tions for positive z are mode 2 failures. Hence slope and curvature are difficult 
to define. The naive normal theory approximation ignores this. It is not clear 
why the asymptotic estimate should tend to be small. 

These results indicate that the approach to asymptotic theoi^ is rapid if the 
model is "estimable" but slow otherwise. They also indicate that rapid approach 
to normal theory approximation is quite possible. Much further simulation is 
required to substantiate these conjectures. 
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7.0 CONCLUSIONS AND HI.' COMMENDATIONS 

Previous studies have concentrated on the competing failure mode estimation 
problem where many samples have been obtained from one single population or at 
most where the location parameter depends on one stress variable and repeated 
observations are made at each level of the stress variable. This paper is oriented 
toward generalizing the multiple linear regression situation where more than one 
stress variable exists and not many replications are performed at each combination 
of stress variables. 

This study has indicated the maximum likelihood approach to be feasible. 

For moderate sample sizes the estimators approach the asymptotic covariance 
structure and asymptotic unbiasedness under certain conditions. When estimating 
parameters for mode i when most of the failures are by mode i the location pa- 
rameters are nearly unbiased while the scale parameter can be substantially biased 
for small sample sizes. Both normal theory and asymptotic theory covariance 
structures appear to hold for small samples. More extensive simulations are re- 
quired to concretely verify and quantify these results. 

In the situation where there is substantial censoring of mode i failures by 
other modes, care must be taken to insure that the proposed model is estimable. 
Otherwise serious biasing and much increased covariances result. 

Based on the results obtained, the following procedure is recommended. 

(1) Perform an ordinary multiple linear regression analysis of all observed 
lifetimes regardless of mode of failure. This will provide rough estimates of the 
terms needed in the response function and an initial guess for the scale parameters. 

(2) Fit the model obtained from all the data to each mode separately using as 
observations only the failures by that mode. This will indicate whether or not all 
terms are estimable. If there is serious ill- conditioning, some terms will have 
to be dropped to provide a model for which all terms are estimable. The ordinary 
least squares regression analysis of this model will then provide initial estimates 
of the location parameter values and the scale parameter. 

(3) Refine these estimates by obtaining the maximum likelihood estimators. 
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This report has not addressed the questions of hypothesis testing and confi- 
dence limits. If the approach to asymptotic behavior Is Indeed as rapid as Indi- 
cated here, then perhaps tests and confidence limits based on the estimate of 
the Fisher Information matrix provided by die sample will perform reasonably. 
Additional simulation studies to Investigate this should prove valuable. 
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TABLE I. - SUMMARY OF RESULTS OF SIMULATION EXPERIMENT A 


(Result* ar* baaed on 1000 simulation*. Standard devia tion o f mean la 
estimated a* observed standard deviation divided by VlOOO.) 
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TABLE n. - MEANS AND ESTIMATED 8TANDABD DEVIATIONS 
OF MEANS OF PARAMETER ESTIMATES IK SIMULATION 
EXPERIMENT B 


(Reoulta are baaed on 100 almulatlooa (or aach cut.] 
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TABLE m. - OBSERVED, ASYMPTOTIC, AND NORMAL THEORY COVARIANCES 
OBTAINED FROM EXPERIMENT B 
[Result* art b***d on 100 almulstlon* of Mch caso.) 






































































TABLE IV. . SIMULATION RESULTS FOR EXPERIMENT C 


|Rmu 1U art baaod upon 100 simulation* of aach cam.) 



































TABU V. - OBSERVED, ASYMPTOTIC, AND NORMAL THEORY COVARIANCES 

FOR EXPERIMENT C 


(•MB - 1) 


[RmuIU b«Md oo 100 oiBulttiooi of oaeh caac.) 
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TABLE V. - Continued, 
(b) <n«2) 
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TABLE V. - Continued. 


(c) (n » 4) 
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TABLE V. - Concluded, 


<d> (n - 8) 















NORMAL AND EXTREME VALUE DENSITY 
FUNCTIONS! 
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Figurt 2. - Tht imms of m ind 8 is function of t for tot ' 
th r»* simpli slits n • 5. 10 20. Sirs I ndiatt on# standard 
trrorolthtmtin. 
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Figure 31 - The location parameters of modes one 
and two as a function of Zf. The dashed lines 
indicate the upper and lower quartiles of each 
distribution. 
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Figured. - The function y^iz) for models A-L 
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Figure 5. - Contlruiti 
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